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ABSTRACT 

We extend our previous studies aimed at modeling Lyman Alpha Emitters (LAEs) to 
the second population of high redshift sources, Lyman Break Galaxies (LBGs), with 
the final goal of investigating the physical relationship between them. We use a set 
of large (w lO^Mpc"^) cosmological SPH simulations that include a detailed treatment 
of star formation, feedback, metal enrichment and supernova dust production; these 
same simulations have already been shown to successfully reproduce a large number 
of observed properties of LAEs (Dayal et al. 2010). We identify LBGs as galaxies 
with an absolute ultraviolet (UV) magnitude Muv ^ consistent with current 

observational criteria. We then compute the evolution of their (a) UV Luminosity 
Function (LF), (b) UV spectral slope, /?, (c) stellar mass function, (d) star formation 
rate (SFR) density, (e) (specific) Star Formation Rate (sSFR), and compare them 
with available data in the redshift range 6 < 2: < 8. With no further parameter tun- 
ing, the model reproduces the redshift evolution of the LBG UV LF, stellar mass 
function, and SFR density extremely well. It predicts a z-independent w —2.2, 
in agreement with the most recently updated data sets at 2; w 6, but in slight ten- 
sion with the steeper /? — Mjjv observed at z « 7. The mean LBG sSFR increases 
from 6.7 Gyr~^ at z 6 to 13.9 Gyr~^ at z w 8, and it is largest for the smallest 
{M^, < lO^'^M©) LBGs, consistent with the recent findings of Schaerer & de Barros 
(2010), and possibly resulting from a downsizing process. From a comparison of the 
simulated LAE and LBG populations, we find no appreciable differences between them 
in terms of the mass- weighted stellar masses/ages, SFR, and dust content; only the 
faintest LBGs with Mjjv ^ — 18(— 19) z « 6(8) do not show an observable Lyman 
Alpha (Lya) line. LAEs hence represent a luminous LBG subset, whose relative ex- 
tent depends only on the adopted selection criteria. For example, using the Lyman 
Alpha equivalent width selection threshold EW > 55A of Stark et al. (2010), the 
LAE fraction increases towards the faintest LBGs. However, for the canonical value of 
EW > 20A , all LBGs with -23 < Muv < -19 would be identified as LAEs at z « 6; 
the fraction of LBGs showing a Lya line decreases with increasing redshift from z « 6 
to 8 due to the combined effects of dust and reionization. We conclude with a brief 
critical model discussion, which emphasizes the uncertainties inherent to theoretical 
EW determinations. 
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1 INTRODUCTION 

We are in a golden age for the search for high redshift 
galaxies, located at the beginning of the cosmic dawn. This 
progress has been made possible by a combination of state of 
the art instruments and sophisticated selection techniques. 
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The sensitivity of instruments such as the recently installed 
Wide Field Camera 3 (WFC3) on the Hubble Space tele- 
scope (HST), the Subaru and Keck telescopes, have been 
instrumental in pushing to faint magnitudes, in an effort to 
search for the first galaxies. This has been complemented by 
refined selection methods including the dropout-technique 
introduced by Steidel et al. (1996); the band in which a 
galaxy 'drops-out' of visibility is used as an indicator of 
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the Lyman Break at low redshifts and the break bluewards 
of the Lya line at high-redshifts (corresponding to 912 A 
and 1216 A in the rest frame respectively), to estimate the 
galaxy redshift. Though this method has been extremely 
useful in identifying high-z LBGs, it has the drawback that 
the exact source redshift cannot be estimated with complete 
confidence. An alternate strategy for finding high-z galaxies 
is the narrow-band technique (e.g. Malhotra et al. 2005; Shi- 
masaku et al. 2006; Kashikawa et al. 2006; Ifu et al. 2010) 
that is based on looking for the relatively unambiguous Lya 
emission, and is often combined with a drop-out criterion us- 
ing broad-band colors (see eg. Ouchi et al. 2008, 2010). The 
foremost advantage in using such galaxies, called LAEs, is 
that the observed wavelength of the Lya line can be used to 
place stringent constraints on the source redshift. 

We start by summarizing the data collected for LBGs 
that formed during the first billion years after the Big Bang, 
i.e. z ^ 6. Using NICMOS data from the UDF, ACS and 
GOODS fields, Bouwens et al. (2007) have identified 627 i- 
dropouts, corresponding to z ~ 6, down to Muv ~ —17. 
Bouwens et al. (2008) have identified 8 z ~ 7 2-dropout can- 
didates, and have found no z « 9 J-dropouts. Using the 
recently installed WFC3 on the HST, Oesch et al. (2010) 
have detected 16 z-dropouts down to Muv ~ —18 between 
z = 6.5 - 7.5. Bouwens et al. (2010a) have used WFC3 
data to push the observations to z = 8.0 — 8.5 by detect- 
ing 5 y-dropouts. Bouwens et al. (2011a) have also found a 
z ~ 10 J-dropout; if confirmed, this would be the farthest 
galaxy known as of date. Using the same WFC3 data but in- 
stead, performing a full SED fitting to the optical-f infrared 
photometry of all the high-z galaxy candidates detected at 
> 5 — (7 in at least one of the WFC3/IR broad-band filters, 
McLure et al. (2010) find a much larger number of candi- 
dates, as expected. However, they note that about 75% of 
the candidates at z > 6.3 and 100% of the candidates at 
z > 7.5 allow for a z < 2 interloper solution. Analyzing the 
three HUDF and deep wide area WFC3 early release data, 
Bouwens et al. (2010b) have found 66 candidates at z = 7 
and 47 candidates at z = 8. Using data from the HAWK-I 
on the VLT, Castellano et al. (2010) have detected 15 z ~ 7 
candidates. Finally, again using the HAWK-I, Laporte et 
al. (2011) have obtained a sample of 10 galaxy candidates 
at z ^ 7.5, lensed by the cluster A2667; of these, they ex- 
pect at least one, and upto 3 galaxies to be genuine high-z 
sources. The most solid piece of information that has been 
collected by such observations is the number density of LBGs 
as a function of Muv, i.e. the LBG UV LF. It is a notable 
success that the UV LFs collected by these different groups, 
using various instruments and selection techniques, converge 
extremely well for any given redshift. 

In spite of this strong observational push, surprisingly 
little effort has been devoted to theoretically modeling 
high-z galaxies. Numerical simulations attempting to model 
galaxy populations and their evolution at z ^ 6 are scarce, 
with the exception of the works by Nagamine et al. (2006), 
Finlator, Dave & Oppenheimer (2007), Zheng et al. (2010, 
2011) and Forero- Romero et al. (2011). Recently, Salvaterra, 
Ferrara & Dayal (2011) have used high resolution simula- 
tions implemented with a careful treatment of metal enrich- 



ment and supernova (SN) feedback, to get hints on the faint 
end slope of the LBG LF {Muv ^ -19.5) between z ^ 5-W 
and put constraints on the physical properties of such faint 
LBGs. The resolution needed to study the faint end, how- 
ever, meant that only a relatively small cosmological volume 
of (10/i~^comoving Mpc; cMpc)'^ could be simulated. 

In addition to the nature of these pristine galaxies, an 
urgent question that needs to be addressed is the relation- 
ship between LBGs and LAEs, which remains ambiguous as 
a result of the different observational selection techniques 
used to detect these two high-z galaxy populations. Gawiser 
et al. (2006) suggest LAEs are less massive and less dusty 
compared to LBGs, and that LAEs represent the beginning 
of an evolutionary sequence in which galaxies increase their 
mass through mergers and dust through star formation. This 
is in agreement with the work of Pentericci et al. (2007, 
2010): although they find LAEs to be less massive and less 
dust enriched with respect to LBGs, they find no appreciable 
difference in the age of LBGs with/without Lya emission. 
However, this is in contrast with the results found by Kornei 
et al. (2010); these authors find LAEs to be older, less dusty 
and lower in star formation compared to LBGs. They con- 
clude that LAEs represent a later stage of galax;y evolution 
in which SN-driven outflows have have reduced the inter- 
stellar medium (ISM) dust covering fraction. Such simplified 
pictures are complicated by results found by other groups 
which are now summarized: Finkelstein et al. (2009) find a 
range of dust extinction, between E{B — V) = 0.5 — 4.5 for 
14 LAEs at z^ 4.5. Nilsson et al. (2009) find that LAEs oc- 
cupy a wide range in color space, implying that not all LAEs 
are young, pristine objects. Lai et al. (2008) have collected a 
sample of 70 LAEs at z 3.1, about 30% of which have been 
detected in the Spitzer 3.6^m band. These IRAC (InfraRed 
Array Camera) detected LAEs are older and more massive 
than the IRAC undetected sample, which leads these au- 
thors to suggest that the IRAC detected LAEs may be a low 
mass extension of the LBG population. Further, Ouchi et al. 
(2008) have discussed the overlap of the LAE and LBG UV 
LFs, pointing out that the overlap between the LFs depends 
on the the EW criterion used to select these two populations. 
Along the same lines, Verhamme et al. (2008) propose that 
at a given redshift, LAEs and LBGs are the same population 
above a given limiting magnitude; for galaxies fainter than 
such magnitude, LAEs are less massive compared to LBGs. 
Finally, Dijkstra & Wyithe (2011) point out that observa- 
tionally, LAEs are chosen using sophisticated color-color se- 
lection techniques; they caution that simple cuts in EW and 
UV magnitude can lead to uncertainties in the simulated 
LAE number densities. 

Theoretically too, the picture has remained equally un- 
clear. While considerable efforts have been devoted to using 
z ^ 5 LAEs as probes of reionization and galaxy evolution 
(e.g. Santos 2004; Dijkstra et al. 2007ab; Kobayashi et al. 
2007; Zheng et al. 2010, 2011; Dayal et al. 2008, 2009, 2010, 
2011; Dayal & Ferrara 2011), only scant effort has been spent 
in establishing a physically motivated connection between 
LAEs and LBGs at these early epochs. One such effort has 
been presented in Shimizu & Umemura (2010) who propose 
that two kinds of galaxies could be visible as LAEs: early 
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starbursts in young galaxies, or delayed starbursts in evolved 
galaxies. These authors further claim that LBGs showing a 
strong Lya line could possibly include the latter type of 
LAEs. 

In this work, our first aim is to simulate the galaxies 
that populate the entire observed LBG UV LF at z ~ 6, 7, 8, 
ranging from galaxies as faint as Muv ~ — 18 to those as 
bright as Muv ~ —23, to reproduce their observed data 
sets including the UV LFs, /3, the stellar mass function, the 
sSFR and the SFR density, thereby pinning down their elu- 
sive physical properties. To this aim, we use state of the 
art cosmological simulations with the same treatment of the 
metal enrichment and SN feedback, as in Salvaterra, Ferrara 
& Dayal (2011). However, to capture the large luminosity 
range required, we use much larger simulation boxes, such 
that the volume probed is (75/i~^cMpc)'^. 

Secondly, the simulations used in this work have already 
been used extensively (Dayal et al. 2009; Dayal, Ferrara & 
Saro 2010) to study the nature of z ^ 6 LAEs and have 
successfully reproduced a number of LAE data sets includ- 
ing the Lya and UV LFs, their spectral energy distributions 
(SEDs) and Lya EWs, to name a few. In this work, our 
aim is to implement this same physically motivated, self- 
consistent model for all the simulated galaxies at z ~ 6, 7, 8, 
to identify the galaxies that would be observationally cho- 
sen as LAEs and LBGs using the narrow-band and drop- 
out techniques respectively, at these redshifts. Then, for the 
first time, we would be able to compare the physical proper- 
ties of LBGs and LAEs, identified from the same underlying 
galaxy population, to study the connection between these 
two galaxy classes in a self-consistent framework. 



2 MODEL 

The simulations used in this work have been carried out us- 
ing a TreePM-SPH code GADGET-2 (Springel 2005) with 
the implementation of chemodynamics as described in Tor- 
natore et al. (2007); interested readers are referred to Tor- 
natore et al. (2010) for complete details. The adopted cos- 
mological model corresponds to the ACDM Universe with 
flm = 0.26, = 0.74, fib = 0.0413, = 0.95, Ho = 73km 
s~^ Mpc~^ and as — 0.8, thus consistent with the 5-year 
analysis of the WMAP data (Komatsu et al. 2009). The 
simulation has a periodic box size of 75h~^ cMpc and starts 
with 512^ dark matter (DM) and gas particles each; the 
masses of the DM and gas particles are ~ 1.7 X 10® h'^MQ 
and ~ 4.1 X lO'^/i^^Mo, respectively. The run assumes a 
metallicity-dependent radiative cooling (Sutherland & Do- 
pita 1993) and a uniform z-dependent UV background pro- 
duced by quasars and galaxies as given by Haardt & Madau 
(1996). The code also includes an effective model to describe 
self-regulated star formation in a multi-phase ISM and a 
prescription for galactic winds triggered by SN explosions, 
Springel & Hernquist (2003). In their model, star forma- 
tion occurs due to collapse of condensed clouds embedded 
in an ambient hot gas. Stars with mass larger then 8Mq 
explode as supernovae and inject energy back into the ISM. 
The relative number of stars of different mass is computed 



for this simulation by assuming a standard Salpeter initial 
mass function (IMF) between 1 and 100 Mq. Metals are pro- 
duced by SNII, SNIa and intermediate and low-mass stars in 
the asymptotic giant branch (AGB). Metals and energy are 
released by stars of different masses by properly accounting 
for mass-dependent lifetimes as proposed by Padovani & 
Matteucci (1993). The metallicity-dependent stellar yields 
have been taken from from Woosley & Weaver (1995) and 
the yields for AGB and SNIa from van den Hoek & Groe- 
newegen (1997). Galaxies are recognized as gravitationally 
bound groups of star particles by running a standard friends- 
of-friends (FOF) algorithm, decomposing each FOF group 
into a set of dis-jointed substructures and identifying these 
using the SUBFIND algorithm (Springel et al. 2001). After 
performing a gravitational unbinding procedure, only sub- 
halos with at least 20 bound particles are considered to be 
genuine structures, Saro et al. (2006). For each "bona-fide" 
galaxy in the simulation snapshots at 2 « 6, 7, 8, we compute 
the mass- weighted age (t*), the total halo/stellar/gas mass 
(Mh/Mt/Mg), the SFR (M.), the mass weighted gas/stellar 
metallicity (Zg/Z^), the mass- weighted gas temperature and 
the half mass radius of the dark matter halo. 

We start by recapitulating the theoretical model that 
was used to identify LAEs using the same simulation snap- 
shots used in this work and the interested reader is referred 
to Dayal et al. (2008, 2009), Dayal, Ferrara & Saro (2010) 
and Dayal & Ferrara (2011) for complete details. The SED 
for each galaxy in each snapshot is computed using the pop- 
ulation synthesis code STARBURST99 (Leitherer et al. 1999), 
using the values of M«, Z, and t, of the galaxy under con- 
sideration; this includes the contribution both from stellar 
and nebular emission. We assume a neutral hydrogen (H i ) 
ionizing photon escape fraction of fesc = 0.02 (Gnedin et 
al. 2008) and hydrogen case-B recombination to calculate 
the intrinsic Lya luminosity, L™'. The intrinsic continuum 
luminosity, is calculated at rest-frame wavelengths 

A = 1350, 1500, 1700 A at z « 6, 7, 8 respectively; though 
these wavelengths have been chosen for consistency with ob- 
servations, using slightly different values would not affect the 
results in any sensible way, given the flatness of the intrinsic 
spectrum in this short wavelength range. For each galaxy, 
the intrinsic Lya EW is calculated as EW'"' = L^"7Lj,"*. 

Both the Lya and continuum photons so produced are 
susceptible to absorption by dust grains present in the ISM; 
only a fraction (/c) of the Lya (continuum) photons pro- 
duced emerge out of the galaxy undamped by dust. Dust 
is produced both by SN and evolved stars in a galaxy. 
However, several authors (Todini & Ferrara 2001; Dwek et 
al. 2007, Valiante et al. 2009) have shown that the contri- 
bution of AGB stars becomes progressively less important 
towards higher redshifts {z ^ 5.7) since the typical evolu- 
tionary time-scale of these stars (^1 Gyr) becomes longer 
than the age of the Universe above that redshift. We there- 
fore make the hypothesis that the dust present in galax- 
ies at z ^ 6 is produced solely by SNII. The total dust 
mass present in each galaxy is then computed assuming: 
(i) 0.5 Mq of dust is produced per SNII (Todini & Ferrara 
2001; Nozawa et al. 2007; Bianchi & Schneider 2007), (ii) 
SNII destroy dust in forward shocks with an efficiency of 
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Table 1. Average physical properties of LBGs. For each of the redshifts considered in this work (col 1), we show the average stellar mass 
(col 2), the average SFR (col 3), the average age (col 4), the average stellar metallicity (col 5), the average intrinsic UV spectral slope, 
i.e. without considering dust attenuation (col 6), the average UV spectral slope including dust attenuation (col 7), the specific SFR (col 
8) and the average color excess (col 9). 
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about 40% (McKee 1998; Seab & ShuU 1983), (iii) a homo- 
geneous mixture of gas and dust is assimilated into further 
star formation, and (iv) a homogeneous mixture of gas and 
dust is lost in SNII powered outflows. 

To transform the total dust mass into an optical depth 
to UV continuum photons, r, we use the results obtained 
for LAEs: for each galaxy, we assume that dust is made up 
of carbonaceous grains and spatially distributed in a slab 
of radius rd = (0.6, 1.0)r, at z ~ (6,7), respectively. Here, 
r« is the radius of the stellar distribution calculated by ex- 
trapolating the results of Bolton et al. (2008) who have de- 
rived fitting formulae to their observations of massive, early 
type galaxies between z = 0.06 — 0.36. The continuum lu- 
minosity that reaches the observer, Lc, is then related to 
the intrinsic luminosity as Lc = Ll"'*fc. In order to achieve 
the best fit to the LAE Lya LP, we slightly modify the 
relation between fa and fc presented in Dayal, Ferrara & 
Saro (2010): for each simulated galaxy at z ~ (6, 7) we use 
U = (1.5,0.6)/ce-'^^'"'*^'= and Mk is independent of red- 
shift with a value lO^^'^M©; we assume the same ratio of 
/a //c at 2 ~ 8 as at 2 ~ 7. Interested readers are referred 
to Dayal, Ferrara & Saro (2010) for complete details of this 
calculation. 

Unlike continuum photons, the Lya photons that es- 
cape out of the galactic environment are further attenu- 
ated by the H i in the intergalactic medium (IGM), and 
only a fraction < < 1 reaches the observer. In turn 
the H I number density at a given redshift depends on the 
assumed reionization history, which we have taken accord- 
ing to the Early Reionization Model (ERM; Gallerani et al. 
2008a). Such reionization history explains a number of LAE 
and QSO observations (see Dayal et al. 2008; Gallerani et 
al. 2008ab), which are used to fix the H i fraction xhi ~ 
(6 X 10"^2.4 X 10"^0.22) a.t z f=i (6,7,8) respectively. We 
then compute the radius of the H ii region each galaxy is 
able to ionize around itself depending on its age and SFR, 
using fesc ~ 0.02 as mentioned before. However, in reality, 
multiple galaxies generally contribute ionizing photons to 
the same ionized region due to source clustering. The size of 
this "collective" H ii region is therefore larger and the pho- 
toionization rate inside it includes the direct radiation from 
the galaxy under consideration, those clustered around it 
and residing in the collective HII region, and the UV back- 
ground. This procedure, described in detail in Dayal & Fer- 
rara (2011), was carried out for each galaxy in the simu- 



lated volume. Assuming photoionization equilibrium within 
the effective ionized region of each galaxy and forcing xhi 
to attain the assigned global value at the edge of this region, 
we used the Voigt profile to calculate the optical depth, and 
hence Ta for Lya photons along the line of sight. The ob- 
served LyQ luminosity for each galaxy is then calculated as 

r int r rp 

2.1 Identifying LAEs and LBGs 

When compared to the standard Sheth-Tormen mass func- 
tion (Sheth & Tormen 1999), the simulated mass function 
is complete for halo mass Mh ^ 10^" M0 in the entire red- 
shift range of interest, i.e. 2 ~ 6 — 8; galaxies above this 
mass cut-off are referred to as the "complete sample". At 
each simulated redshift, z ^ 6,7, 8, in consistency with the 
current observational criteria, galaxies from the "complete 
sample" with an absolute magnitude Muv ^ — 17 are iden- 
tified as LBGs; galaxies with La ^ lO^^ '^erg s"'^ and an 
observed Lya equivalent width EW — La/Lc ^ 20 A are 
identified as LAEs. 



3 COMPARISON WITH LBG OBSERVATIONS 

We reiterate that galaxies identified as LAEs in this work are 
the same as those presented in Dayal et al. (2009) and Dayal, 
Ferrara & Saro (2010). These have already been shown to 
successfully reproduce a number of LAE data sets including 
the Lya/continuum LFs, EW distributions, observed SEDs 
and the color excess, E{B — V), to mention a few. We re- 
fer the interested reader to the above mentioned works for 
complete details. In this work, we limit ourselves to com- 
paring the results for the theoretical LBGs to observations. 
We start by comparing the theoretical LBG UV LFs to the 
observed ones, for each of the redshifts of interest. 

3.1 UV Luminosity Functions 

As mentioned before, the LBG UV LFs are by far the most 
solid piece of information collected for high-z galaxies. In- 
deed, the observations by different groups are all in excellent 
agreement at z ~ 7, 8, as seen from Panels (b) and (c) of 
Fig. [T] Reproducing these LFs is therefore the first test of 
our model: for each redshift of interest, the theoretical LF 
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Figure 1. LBG UV LFs for z fii 6,7, 8 as marked in eacii panel. In all panels, points sliow the observed data, solid lines show the 
theoretical results and the shaded regions show the poissonian errors associated with the latter. The observed LBG UV LFs have been 
taken from: (a) z 6: Bouwens et al. (2007; filled circles); (b) z ^ 7: Oesch et al. (2010; filled squares), Bouwens et al. (2010a; empty 
circles), Bouwens et al. (2010b; filled circles), Castellano et al. (2010; empty triangles) and McLure et al. (2010; filled triangles); (c) 
z 8: Bouwens et al. (2010a; empty circles), Bouwens et al. (2010b; filled circles) and McLure et al. (2010; filled triangles). 



is computed by counting the number of LBGs in each Muv 
bin and dividing it by the simulated volume, and the size of 
the bin. 

We find that at all the redshifts considered, both the 
amplitude and the slope of the theoretical LFs are in excel- 
lent agreement with the observations, over a broad magni- 
tude range Muv ~ — 18 to —23, as seen from panels (a,b,c) 
of Fig. [1] We consider this to be an encouraging success of 
our model, particularly because once L*"* and fc are cal- 
culated using the physical properties of each galaxy, as ex- 
plained in Sec. (2] there are no free parameters left to calcu- 
late the UV LF. 

We briefly digress here to discuss the LAE UV LFs: we 
begin by reiterating that, as explained in Sec. [21 the model 
parameters used for both LAEs and LBGs have the same 
values at a given redshift. The parameters used in this work 
have already been shown to reproduce the LAE UV LFs ex- 
tremely well at z ~ 6, 7 (see Fig. 7 of Dayal, Ferrara & Saro 
2010). Indeed this is not a surprise since observationally, 
both the amplitude and shape of the LAE UV LFs coincides 
perfectly with that observed for LBGs between Muv ~ —22 
to —20, at 2 « 6, 7; this agreement has already been pointed 
out and discussed in detail in Shimasaku et al. (2006). Such 
agreement hints at the fact that LAEs and LBGs are simi- 
lar galaxy populations (at least in the above UV magnitude 
range) . 

We now return to the discussion regarding the LBG 
LFs where we find that the theoretical LFs shift towards 
lower luminosity or decreasing Muy with increasing red- 
shift, mimicking a pure luminosity or density evolution. 
The value of the best-fit Schechter parameters to the the- 
oretical LFs can be quantified by a characteristic mag- 
nitude Muv ^ (-20.3 ± 0.1,-20.1 ± 0.2,-19.85 ± 0.15) 
and slope a ^ (-1.7 ± 0.05,-1.7 ± 0.1,-1.6 ± 0.2) for 
z « (6, 7, 8) which, within errors, are consistent with the 
values Muv ~ (-20.24 ± 0.19, -20.14 ± 0.26, -20.1 ± 0.52) 



and Q ^ (-1.74 ± 0.16, -2.01 ± 0.21, -1.91 ± 0.32) derived 
by Bouwens et al. (2010b) for the same redshifts. 

Note, however, that the large cosmological volume simu- 
lated (« 10® cMpc^) naturally limits the mass resolution. As 
mentioned in Sec. [51 when compared to the standard Sheth- 
Tormen mass function, the theoretical halo mass function 
is only complete for galaxies with halo mass ^ 10^" M© 
at any of the redshifts considered here. We are thus un- 
able to resolve the sources populating the faint end of the 
LF, i.e. Muv ^ —18, which explains the flattening of all 
the theoretical LFs below this magnitude limit, as shown in 
Fig. [T] As a consequence, we are also unable to constrain 
the faint end slope of these LFs. Readers interested in the 
faint end of the LF are referred to a previous work (Sal- 
vaterra, Ferrara & Dayal 2011) where the authors have used 
higher resolution simulations with an identical implementa- 
tion of physical processes, albeit with a much smaller volume 
(lO'^cMpc''), to flnd an almost constant value of a ~ —2 for 
the faint end slope between 2 5 — 10. 



3.2 UV spectral slopes 

The UV spectral slope, /?, of galaxies is generally parame- 
terized using a power law index such that the specific flux 
fx oc A'' (e.g. Meurer et al. 1999); in this convention, /? — —2 
corresponds to a source that has a zero color in the AB mag- 
nitude system. Since P depends on the intrinsic properties 
of the galaxy, it has become popular to get a hint on the 
values of t,, Z. and E{B - V) for high-z LBGs. 

However, the values of /? inferred from broad band pho- 
tometry and their evolution with redshift, have been hotly 
debated over the past few months as is now summarized. 
Using HUDF and ERS data, Bouwens et al. (2010a) find 
smaller /? (i.e. bluer spectra) both with decreasing magni- 
tude and increasing redshift: LBGs at 2 ~ 7 are in general 
bluer than those at z ~ 4, and for 2 ~ 7 LBGs, /3 de- 
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Figure 2. The LBG UV spectral slope, /3, as a function of the UV magnitude (in bins of 1.0 dex) at z ^ 6, 7, 8, as marked in each panel. 
Empty (filled) squares show the theoretical /3 values without (with) dust correction with the error bars showing the 1 — cr error. Filled 
circles show the /3 values inferred by Dunlop et al. (2011) for their 'robust' sample, with the additional condition of at least one 8 — cr 
near-IR detection (see Tab. 1 of Dunlop ct al. 2011), empty circles show the results from Bouwens et al. (2011b), with the associated 
random and systematic errors and filled triangles show the f) values inferred by Finkelstein et al. (2011). 



creases from ~ —2 to ~ —3 as Muv increases from —20.5 
to —18.5. Dunlop et al. (2011) have used the same data 
set to explore the efTects of including/excluding less robust 
LBG candidates and have considered the fact that the fields 
used have different limiting magnitudes, to study the ef- 
fects of noise and selection bias at any given UV luminosity. 
These authors claim that artificially low values of /3 can 
be found in the deepest 0.5 magnitude bin of the WFC3 
selected samples, irrespective of luminosity or redshift (see 
also Finkelstein et al. 2010). Confining their analysis to ro- 
bust LBG candidates, these authors find no trend of P with 
either magnitude or luminosity: their average P value is con- 
sistent with {/?) « —2.05 ±0.1 over z « 5 — 7 and between 
Muv ~ —22 to —18. The results obtained by Dunlop et al. 
(2011) are supported by the recent work of Wilkins et al. 
(2011) who find a zero mean UV continuum color in the AB 
magnitude system, i.e. /3 « —2 between Muv ~ —21 to —18 
for WFCS-selected LBGs at z « 7. Finkelstein et al. (2011) 
have fit spectral synthesis models to the LBG data collected 
from CANDELS, the HUDF and ERS, to infer the best-fit 
value of p. These authors also find no trend of /3 with the 
UV magnitude at 1500A in the rest-frame; their average /? 
values of (fi) ^ {-2.07toM, -2.37l:o^ol) at z ^ (6,7) are 
consistent, within errors, with the results found by Dunlop 
et al. (2011). Recently, Bouwens et al. (2011b) have updated 
their results using data from HUDF09 and CANDELS data, 
to find {/3) « —2.2 at z « 6, 7. These new values are now 
consistent with the results found by Dunlop et al. (2011, 
Wilkins et al. (2011) and Finkelstein et al. (2011) mentioned 
above. 

We start by obtaining the intrinsic value of /3 {Pint; i.e. 
without considering the effects of dust) for each of the sim- 
ulated LBGs on the UV LFs at jz ~ 6, 7, 8, shown in Fig. 
[1] We use 30 values of P evenly spaced between —3 and 
to fit a line through the intrinsic SED of each LBG (as 
obtained using STARBURST99) between 1500 - SOOOA in the 



rest frame, sampled every 100 A; the value of P yielding the 
smallest x^-error is then chosen as the P value for that par- 
ticular galaxy. We find that Pi„t decreases slightly for fainter 
galaxies at any given redshift as shown in Fig. [2] For exam- 
ple. Pint decreases from —2.3 to ~ —2.5 as Muv increases 
from —23.5 to —17.5 at z « 6; this is due to the fact that 
smaller (or less luminous) LBGs are in general younger and 
have lower stellar metallicity, as will be discussed in more 
detail in Sec. [4] below. Moreover, as shown in Tab. [1] on 
average, the values of and of LBGs shift to progres- 
sively lower values with increasing redshift: the average age 
and metallicity shift from U f» 142 Myr {Z, a; O.I4Z0) at 
z « 6 to « 69 Myr (Z. ^ O.OSZq) at z f» 8. All this 
results in a progressively bluer SED with redshift, such that 
{P^nt} ~ (-2.45, -2.5, -2.55) at z ^ (6, 7, 8). 

We then calculate the observed P values, taking into ac- 
count the dust enrichment of LBGs with redshift. We con- 
volve the intrinsic spectrum of each LBG with the SN extinc- 
tion curve (Bianchi & Schneider 2007); the latter is scaled for 
the fc calculated for each LBG (see Sec. [21). We use the SN 
curve since it has been shown to successfully interpret the 
observed properties of the most distant quasars (Maiolino 
et al. 2006) and gamma-ray bursts (Stratta et al. 2007). 
We again use the same fitting procedure described before to 
find the P value minimizing the x^-error. As expected, in- 
cluding dust attenuation makes P redder (see Fig. [21), such 
that the average values are: (P) « (—2.2, —2.28, —2.32) for 
z « (6,7,8). As seen from these values, the difference in 
P with/without dust effects decreases with increasing red- 
shift; this is easily explicable by the fact that galaxies are in 
general younger at increasing redshifts, and therefore have 
less time to produce dust. This argument is supported by 
the decreasing stellar metallicity values shown in Tab.[T] Fi- 
nally, we find that P decreases slightly with increasing mag- 
nitude at any redshift; for example, it changes from ~ —2.1 
to ~ —2.25 as Muv increases from —23.5 to —17.5 at z ~ 6. 
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Figure 3. The LBG stellar mass function at z « 6, 7, 8 as marked in each panel. In each panel, the solid line shows the theoretical result 
(poissonian errors are invisible on this scale); the completeness limit of the simulation is M» Filled (empty) points show the 

corrected (uncorrected) LBG stellar mass function inferred by Gonzalez et al. (2011) with error bars reflecting the uncertainty in the LF 
and the mass-luminosity relation used (see Gonzalez et al. 2011 for details). 



We now compare the above predictions with the most 
recent available data. We start by noting that within error 
bars, both the amplitude and the trend of the theoretical 
P — Muv nicely fits the data sets of Dunlop et al. (2011), 
Bouwens et al. (2011b) and Finkelstein et al. (2011) at z ~ 6. 
At « fti 7, the P — Muv trend inferred observationally by 
these same authors appears to be steeper than the theoreti- 
cal one. Such a steepening of the /3 — Muv relation remains 
hard to understand theoretically, given that the dust enrich- 
ment, ages and metallicities of both LAEs and LBGs do not 
evolve appreciably in z — 6 — 7 (see Sec. [4] that follows). 
It remains to be seen whether improvements in data qual- 
ity/analysis, as it has been the case at 2 ~ 6, can release 
such tension between theory and observations at 2: = 7. 

Finally, as shown in Fig. [2] and discussed above, we 
find that /3 decreases only slightly with increasing red- 
shift/magnitude and is consistent with (/3) ~ —2.2 at 2 « 
6, 7, 8. This result is consistent with all the most updated 
observational results available, including those of Dunlop et 
al. (2011), Wilkins et al. (2011), Bouwens et al. (2011b) and 
Finkelstein et al. (2011). 

3.3 Stellar Mass Functions 

At each redshift, the stellar mass function is calculated by 
counting the number of LBGs in each (logarithmic) stellar 
mass, M« , bin and dividing it by the total simulated volume, 
and the size of the bin. As a first result, we see from Fig. [3] 
that the number density of LBGs in any given mass bin de- 
creases with increasing redshift, as expected in standard hi- 
erarchical structure formation models, where galaxies build 
up their stellar mass gradually from the merger of smaller 
objects. Hierarchical scenarios also naturally explain the 
upper-mass cutoff shifting to progressively lower M, values 
with redshift, as the largest objects, with M* « lO^^'^M© at 
2; « 6, have not had sufficient time to grow at 2 8, where 
the upper mass cut-off is M, f» 1O^'®M0. 



We compare these theoretical stellar mass functions to 
those derived by Gonzalez et al. (2011) by SED fitting of 
the broad band photometric colors obtained from Hubble- 
WFC3/IR observations combined with deep GOODS-S 
Spitzer/IRAC data. This is a good consistency check of our 
model, given that Gonzalez et al. (2011) assume a constant 
star formation, to find the stellar masses that best fit the 
observed SEDs. We find the agreement between the the- 
oretical mass function and the data is quite satisfactory, 
both in terms of the slope and the magnitude of the dust- 
corrected estimates for stellar mass M« ^ IO^'^Mq at any of 
the redshifts considered; varying the model parameters such 
as modelling star formation to occur in bursts, or adding 
the contribution from nebular emission might affect such an 
agreement (see Schaerer & deBarros 2010). Below this mass 
threshold, our simulations start being incomplete because of 
the numerical resolution issues already discussed in Sec. 2; 
the decrease in the mass function is probably an artifact of 
such resolution issues. 



3.4 Star Formation Rate Densities 

Another useful quantity that can be extracted from our nu- 
merical simulations is the LBG SFR density, p.. To com- 
pare directly with available data, at each of the simulated 
redshifts, this value is computed by summing the intrinsic 
SFR of all LBGs whose observed luminosity is larger than 
0.05Z/*=3 (corresponding to Muv,ab = —21.07 at 2 = 3) and 
dividing this sum by the total simulated volume. We start by 
observing that p* decreases with increasing redshift, tracing 
the cosmic SFR density (e.g. see Tab. 2 of Hopkins 2004). 
Further, as seen from Fig. 31 the theoretical SFR density 
values are in good agreement with the dust-corrected (un- 
corrected) estimates of Bouwens et al. (2011b) for 2 « 6 
(2^7,8). 

This agreement, however, is probably the result of an 
IMF-dust attenuation degeneracy as we explain in the fol- 



© 2002 RAS, MNRAS 000, [T]-?? 



Dayal & Ferrara 







O 




-2.5 - 



1.5 



o 1 



- 



T — e — I 



I 



I 1 



Figure 4. The SFR density for LBGs, obtained by integrating 
down the UV continuum luminosity density to 0.05L*_3, plotted 
as a function of redshift. Filled squares show the theoretical re- 
sults with error bars showing the poissonian errors; solid (empty) 
triangles show the dust uncorrected (corrected) values inferred by 
Bouwens et al. (2011b) assuming a Salpeter IMF between 0.1 and 
125 Mq and a constant SFR for > 100 Myr. 



lowing. Bouwens et al. (2011b) have used a Salpeter IMF 
between 0.1-125 M© and fc « (0.66,0.77,0.77) at 2 « 
(6, 7, 8) respectively, as seen from a comparison of their 
dust corrected/uncorrected points at these redshifts (see 
Tab. 7, Bouwens et al. 2011b). Our simulations instead, use 
a Salpeter IMF between 1 - IOOMq; the IMF lower limit 
of IMq results in an intrinsic luminosity which is about 
2.5 times larger compared to that from the IMF used by 
Bouwens et al. (2010b). As expected, to fit the UV LF am- 
plitude, our fc S3 (0.2, 0.3, 0.36) values at z ~ (6, 7, 8) (cor- 
responding to a mean E{B -V)^ (0.16, 0.12, 0.1), see Tab. 
[T]) are then about 2.5 times smaller than those of Bouwens 
et al. (2010b). Thus, although theoretical dust corrections 
also decrease with increasing redshift, they are not negligi- 
ble even in the most distant LBGs (see also Schaerer & de 
Barros 2010 and Laporte et al. 2011). 

We now turn to our predictions for the specific SFR, 
where for each galaxy, the sSFR is calculated as sSFR = 
M«/M«. At each of the simulated redshifts, the value of 
sSFR in a given M, bin is the average over all the LBGs 
in that mass bin. First, we find that the sSFR decreases 
with increasing M*; galaxies with M, = lO^'^'^'^M© and 
Mt > IO^'^Mq show less than half of the sSFR of galaxies 
with M, < 1O*'^M0, at any of the redshifts considered here. 
This apparently puzzling result might be the smoking gun 
of downsizing: whereas more massive galaxies formed large 
amounts of stars in a short time interval which quenched 
later star formation, dwarf galaxies built up their stellar 
budget more gradually. Thus, by the time the SFR is roughly 



Figure 5. Specific SFR for LBGs plotted as a function of redshift. 
Filled squares show the theoretical results averaged over all stellar 
masses with the shaded region showing the associated poissonian 
errors; dotted, short and long dashed lines show the theoreti- 
cal estimates for galaxies with total stellar mass Af» < 
M, = 10*'5-9-5Mq and M, > respectively. Empty sym- 

bols refer to observations: the star shows the estimate of Yabe et 
al. (2009) for galaxies with mass between 10*~^^Mq with a me- 
dian mass of 4.1 X IO^Mq. Circles show the estimates from Stark 
et al. (2009). The triangle shows the average value from Gonza- 
lez et al. (2010) whose sample has a median mass of 5 X IO^^Mq. 
Finally, the upper and lower empty squares show the upper and 
lower limits inferred by Schaerer & de Barros (2010) and refer to 
galaxies with masses ^ 10* and 10^' ^Mq respectively. 



similar for the two populations, more massive galaxies al- 
ready contain a large amount of stellar mass, and are there- 
fore characterized by a lower sSFR. This trend appears to 
hold at all redshifts, as seen from Fig.[5l We note that our es- 
timates of the sSFR of galaxies with M* = 10*'^"^-^Mq are 
consistent with the those obtained by Schaerer & de Barros 
(2010) for galaxies with M. = lO^ '^M© (Fig. [5J. Further, 
averaged over all LBGs at a given redshift, the sSFR rises 
with redshift, from « 6.7 to ~ 13.9 Gyr"'^ from z « 6 to 
z « 8, as shown in Fig. [5] and quantified in Tab. [T] Such 
a trend is easily explained combining the fact that galaxies 
shift to progressively lower halo (and stellar) masses with 
redshift (see Fig. |3} and as mentioned above, the sSFR rises 
with decreasing stellar mass. 



4 LAE-LBG CONNECTION 

Having identified the simulated galaxies that would be ob- 
servationally selected as LBGs, and shown that they repro- 
duce a number of observed data sets, we can now compare 
their properties with those of the LAEs previously found in 
the same simulations. 
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Figure 6. As a function of Muv (in bins of 1 dex), we sliow tlie pliysical properties of LBGs (squares) and LAEs (triangles). The 
physical quantities represented by the rows are (a) top row: the total stellar mass, (b) middle row: the mass weighted stellar age and (c) 
bottom row: the total dust content. Columns from left to right refer to z Ri 6, 7, 8 as marked and the error bars show I — a errors. 



We start by considering the stellar mass, Af*, which 
is possibly the best constrained physical quantity available 
for high-z galaxies when near-IR data is available. We find 
that for galaxies brighter than Muv ^ ( — 18, —19, —19) at 
z ~ (6,7,8), LAEs and LBGs show very similar values of 
Mt in a given Muv bin. The stellar mass and UV mag- 
nitude are tightly coupled at each redshift; galaxies with 
the largest M» (and hence M«; see Fig. 7 of Dayal et al. 
2009) are the most UV luminous, implying that although 
dust obscuration can modify the slope of the SFR-UV lu- 
minosity relation, it cannot alter its monotonic trend. A 
persistent evolutionary feature is that the LBG popula- 
tion extends to lower stellar masses/fainter luminosities: 
this is because the smallest galaxies do not produce enough 
Lya to be visible as LAEs, according to the adopted se- 
lection criterion {La ^ lO'^^'^ergs"^ and EW ^ 20 A). 
Finally, the M« ranges shift to progressively lower values 
with redshift (see also Sec. I3.3|l for both LAES and LBGs, 
as expected in a hierarchical structure formation scenario; 
at 2 ~ 6 the stellar masses for LBGs (LAEs) range be- 
tween 1O^'^"^"'^(1O*'^"^°-^)M0, while they shift down to 
M. ^ lO*-^'^(lO**-''-^ *)M0 at z 8 (Fig. |6}. The LAE 
stellar mass range evolves faster than that for LBGs, due to 



the combined effects of dust and reionization, as explained 
in Sec. [5] below. 

Similar considerations hold for the stellar ages, t* , of the 
two populations, as they are found to be very similar in any 
given absolute magnitude bin for AIuv ^ ( — 18, —19, —19) 
at z ~ (6, 7, 8), further pointing to the tight relation between 
these two types of galaxies. Both LBGs and LAEs are com- 
posed of intermediate age systems, with 100 < t*/Myr < 
250 (60 < i*/Myr < 250) at z « 6 (8), as seen from pan- 
els (d)-(f) of Fig. |6] As expected, the average t, decreases 
with redshift; in agreement with the downsizing scenario 
advocated above, larger systems are more dominated by 
older populations than smaller systems. Finally, the galaxies 
brightest in the UV are on average the oldest, by virtue of 
their largest M« values. 

Finally, as shown in panels (g)-(i) of Fig.[Sl both LAEs 
and LBGs have nearly identical dust masses (again, apart 
from the faintest magnitude bins), as a result of their over- 
lapping stellar masses and ages. As galaxies are on av- 
erage younger and less massive at higher redshifts, more 
distant objects tend to be less dusty such that the dust 
mass, Maust = 10'^-''"''(10^'^"^)Mq at z 6 and Mdust = 
;^q5.5-6.5^;^q6.i-6.5^j^^ at z 8 for LBGs (LAEs), respec- 
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Figure 7. As a function of Mi/y (in bins of 1 dex), we show ttie intrinsic (top three panels) and the observed Lya EW (bottom three 
panels). The values for LAEs (triangles) and LBGs (squares) are shown, with error bars referring to 1 — cr errors. Columns refer to results 
for 2 Si 6, 7, 8 as marked. 



tively. Finally, as expected, the dust content of galaxies in- 
creases with M« (decreasing Muv); the brightest galaxies 
have formed the largest amount of stars over a longer pe- 
riod of time compared to smaller systems. 

From the above discussion, we conclude that high- 
redshift sources identified as LAEs and LBGs essentially 
sample the same underlying galaxy population; however, 
LAEs represent a subset of a larger LBG sample in the 
sense that only the faintest LBGs do not show an observ- 
able Lya line. We clarify that LAEs are neither younger/less 
dusty/smaller nor older/more dusty/larger as compared to 
LBGs, but rather a wide range of LBGs show an observ- 
able Lya line. Indirectly, our findings are supported also by 
the recent results of Malhotra et al. (2011) who find similar 
physical sizes for LAEs and LBGs at z ^ 5. Intriguingly, 
studying 92 (z) = 2.65 UV continuum selected galaxies to 
very low surface brightness limits, Steidel et al. (2011) find 
that extended, diffuse Lya emission increases the total Lya 
flux by a factor of about 5 on average; accounting for such 
diffuse emission, essentially all their LBGs would qualify as 
LAEs. 

Finally, we caution the reader that the overlap between 
the LAE and LBG populations depends sensitively on red- 
shift dependent observational selection criterion such as the 
UV magnitude, Lya flux and EW limits, used to select these 
two populations; the results found in this work which sug- 
gest LAEs to be a subset of LBGs might not necessarily hold 
true for data sets at other redshifts selected using different 
selection criterion. This point has already been made in a 
number of previous works: Ouchi et al. (2008) have shown 



that the fraction of LBGs showing a Lya line depends on 
the Lya EW cut imposed at any given redshift, Verhamme 
et al. (2008) have shown that while LAEs brighter than a 
certain limiting magnitude (which is a function of redshift) 
are the same as LBGs, at fainter magnitudes, LAEs are less 
massive compared to LBGs. Also, Dijkstra & Wyithe (2011) 
have pointed out that simulating LAEs using simple cuts in 
EW and UV magnitude lead to uncertainties in the pre- 
dicted LAE number densities. 



5 IMPRINTS OF REIONIZATION 

As explained in Sec. (2] while both Lya and continuum pho- 
tons are attenuated by ISM dust, the IGM ionization state 
only affects the observed Lya luminosity, leaving the ob- 
served continuum luminosity unchanged; a decrease in the 
observed EW with increasing redshift could then possibly 
indicate an increasingly neutral IGM, a test proposed by 
various authors (e.g. Stark et al. 2010). We now discuss both 
the intrinsic and observed Lya EWs for the LAEs and LBGs 
identified at each of the simulated redshifts. 

As discussed in Sec. |4l for a given Muv value, both 
t, and Z^, shift to progressively lower values with increas- 
ing redshift, for both LAEs and LBGs. Since the stellar 
spectra of younger, more metal-poor populations are harder 
(i.e. stars produces more H i ionizing photons) the value of 
EW^"^ increases with redshift for any UV magnitude bin; 
for the same reasons, fainter galaxies exhibit larger EW^"^ 
values at all redshifts. Such trends hold for both LAEs and 
LBGs by virtue of their nearly identical physical proper- 
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Figure 8. As a function of Muv (in bins of 1 dex), we show the fraction of LBGs showing a Lya line. The left and right panels refer to 
the fraction of LBGs showing a Lya line with EW > 55A and EW > 20A , respectively. The points in the left panel show the fractions 
obtained by Stark et al. (2010) in the range 3 < z < 6 and the histogram shows the model results for z ^ (i\ the results for z Ri 7, 8 are not 
visible on this plot. In the right panel, histograms from the widest to the narrowest show the model results for 2 6, 7, 8, respectively. 
See text for details. 



ties, as shown in panels (a)-(c) of Fig. [T] Quantitatively, the 
range of SH^'"' fa 90- 130A (brightest to faintest galaxies) 
at 2 f» 6 and increases to ~ 110 — 180 A at 2 ~ 8. 

Before we discuss the observed EWs, we recall that the 
following relation holds: EW = EW'"\U/ fc)Ta„ where Tc 
is the transmissivity of Lya photons through the IGM. At 
2; « 6 — 7, according to the Early Reionization Model that 
best fits the observations, the H i fraction is xhi ~ 10~^; 
this largely ionized IGM results in transmissivity values con- 
strained in the narrow range Ta ~ 0.4 — 0.6, from the small- 
est to the largest LAEs (see also Fig. 3, Dayal et al. 2009). 
At z ~ 8, instead, the more neutral IGM (xhi ~ 0.22) sub- 
stantially and preferentially damps Lya photons from the 
smallest galaxies which therefore become undetected in nar- 
row band observations; the transmissivity decreases globally 
to « 0.05 - 0.4. 

Further, we remind the reader that have used fa/fc = 
(1.5,0.6,0.6)6"*'^''/*^'= at z « (6,7,8), where Mh is the 
galaxy halo mass and Mfc — IO^^'^Mq, as required to best 
reproduce the amplitudes of the cumulative LAE Lya LFs. 
When taken together, these facts result in the value of 
EW rising smoothly from ^ 30 A at Muv ~ -23.5 to 
EW « 80 A, at Muv ~ -17.5, at z 6. However, the 
value of EW is much fiatter at z ~ 7 as a result of the lower 
escape fraction of Lya photons from the ISM, compared to 
that for the continuum photons (fa/fc ~ 0.6) required to 
match the LAE LFs. Finally, at z ~ 8, the trend reverses, 
with the smallest galaxies showing the smallest EW, as a 
result of their strongly reduced Lya transmission through a 
partially neutral IGM. 

Given the above discussion, is then clear that reioniza- 
tion is the key factor (but by no means the only one) reg- 
ulating the relative fraction of observed LBGs and LAEs. 
With our results at hand, it is straightforward to predict 
the fraction of LBGs that show a visible Lya line (i.e. 
La ^ 10^^'^erg s~^), at any given redshift. An observational 
consensus on this quantity is far from being reached, as it 
will be clear from the following. For example, using Keck 



spectroscopy of about 600 LBGs in the range 3 < z < 6, 
Stark et al. (2010) find that the least luminous LBGs most 
often show Lya emission with EW > 55A. To compare to 
such observations, we select all LBGs that would be iden- 
tified as LAEs with EW > 55A, at each of the redshifts 
considered. We find that at z ~ 6, indeed it is the least lu- 
minous LBGs most often show Lya emission as shown in 
panel (a) of Fig. [S] qualitatively our results broadly agree 
with those of Stark et al. (2010). However, at higher red- 
shifts we find only extremely few (z = 7) or no (z = 8) 
LAEs with EW > 55A due to the low escape fraction of 
Lya and increasing IGM opacity, as explained above. These 
combined effects lead to a negligible fraction of LAEs being 
found amongst z > 6 LBGs in our model. 

However, the most often adopted EW cut corresponds 
to EW > 20A. Using this value, we find that within the 
range Muv ~ —23 to —20, all LBGs would be classified as 
LAEs at z ~ 6. This is perhaps not surprising, given that 
the amplitudes of both the LAE and LBG UV LFs match 
perfectly within this range, as has already been pointed out 
by Shimasaku et al. (2006) and discussed in detail by Ouchi 
et al. (2008) . The magnitude range of LBGs showing a Lya 
line decreases with redshift: the smallest galaxies do not even 
produce enough Lya, as mentioned before; the fraction of 
the largest LBGs showing Lya emission progressively de- 
creases with redshift due to the dependence of fa on Mh ex- 
plained above. Further, the decrease in the fraction of LBGs 
showing Lya emission between z ?a 6 — 8 shown above, has 
recently also been found in observational data by both Pen- 
tericci et al. (2011) and Schenker et al. (2011); these authors 
find that the fraction of LBGs showing Lya emission drops 
between z ~ 6 — 7, which they interpret as a sign of the IGM 
becoming progressively neutral (see also Fontana et al. 2010 
and Vanzella et al. 2011). We reiterate that although we see 
such a trend, between z « 6 — 7 it is driven predominantly 
by the dust distribution (which governs the ratio fa/fc); be- 
tween z ~ 7 — 8 on the other hand, this trend arises due to 
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a combination of the dust distribution, and the IGM being 
more neutral. 

As a final note of caution, we restate that results re- 
garding the Lyof EWs depend sensitively on a number of as- 
sumptions and modeling of poorly understood physical pro- 
cesses, namely: (a) we have used a luminosity-independent 
value of fa/ fc, for all galaxies at a given redshift. This is 
empirically fixed by normalizing the theoretical LAE LPs to 
the observed ones. While a number of works (Neufeld 1991, 
Hansen & Oh 2006) have emphasized that such ratio de- 
pends on the dust distribution (i.e. homogeneous/clumped) 
in galaxies, at the moment, essentially no experimental hint 
is available at high-2 on this aspect; (b) the effects of peculiar 
velocities in calculating Tc, have not been included; infiows 
(outflows) blueshift (redshift) the Lya line, thereby decreas- 
ing (increasing) Ta- A further complication arises due to the 
fact that both these above mentioned effects depend on the 
line of sight between the LAE and the observer. Finally, 
we note that we have used the stellar mass weighted age 
and assumed the SFR to be constant over such time period 
to compute the intrinsic spectrum of each galaxy. However, 
using the complete star formation history of each galaxy 
might affect the results presented here. Each of these points 
is discussed in more detail in the next Section. Alternatively, 
the tension between the model and the observations could 
also arise due to poorly understood biases/selection effects 
present in the data. 



6 SUMMARY AND DISCUSSION 

We have extended our previous studies aimed at modeling 
Lyman Alpha Emitters (LAEs) to the second major pop- 
ulation of high-z sources, Lyman Break Galaxies (LBGs), 
with the final goal of investigating the physical relationship 
between them. In a set of large (~ lO^cMpc"^) cosmologi- 
cal SPH simulations, including a detailed treatment of star 
formation, feedback, metal enrichment and supernova dust 
production, we identify LBGs as galaxies with an absolute 
magnitude Muv ^ —17, in consistency with current obser- 
vational criteria; the same simulations have already been 
shown to match a number of observed data sets for LAEs 
(Dayal et al. 2010). For the LBGs identified at each of the 
redshifts of interest, i.e. z « 6,7,8, we compute the evolu- 
tion of their (a) UV LF, (b) UV spectral slope, /3, (c) stellar 
Mass Function, (d) SFR density, (e) sSFR, and compare 
them with available data. The main results from this study 
are summarized as follows. 

• The shape and amplitude of the observed LBG UV 
LFs at 2: ~ 6, 7, 8 are reproduced extremely well without 
needing to invoke any free parameters in addition to those 
used to model the LAE Lya/UV LFs; the best-fit Schechter 
function values obtained are M^y ~ (—20.3 ± 0.1, —20.1 ± 
0.2, -19.85 ±0.15) and a ^ (-1.7 ±0.05, -1.7 ±0.1, -1.6 ± 
0.2) at 2 ~ 6, 7, 8, respectively. 

• The mean intrinsic LBG UV spectral slope shows lit- 
tle variation, either with magnitude (in AIuv ~ —22.5 to 
-17.5) or redshift, such that (/?'"') = (-2.45,-2.5,-2.55) 
at 2 « (6,7,8). When convolved with the SN extinction 



curve, the observed spectral slope increases such that (/3) = 
(—2.2, —2.28, —2.32) in the same redshift range, again vary- 
ing only slightly between the brightest and faintest galax- 
ies. This is in line with the recent findings of Dunlop et al. 
(2011), Wilkins et al. (2011), Bouwens et al. (2011b) and 
Finkelstein et al. (2011). 

• The LBG stellar mass function matches both the slope 
and amplitude of the values inferred from observations by 
Gonzalez et al. (2010) for M* ^ lO^'^M© at aU of the red- 
shifts considered here; below this mass cut-off, numerical 
resolution limits the accuracy of our prediction. 

• The LBG sSFR rises with redshift such that sSFR^ 
6.75 (13.9) Gyr~^ at 2 « 6 (8), with the smallest galax- 
ies showing the largest sSFR values. We have suggested 
that such behavior could be a manifestation of downsiz- 
ing. Finally, our results agree well with the dust-corrected 
SFR density values inferred by Bouwens et al. (2011b) for 
2 = 6-8. 

• As for the LAE-LBG connection, at 2 « 6, 7, the LAE 
and LBG UV LFs have the same amplitude and slope be- 
tween —22 < AIuv < —20; within this magnitude range all 
LBGs would be identified as LAEs and vice- versa, according 
to the selection criteria usually adopted (see Sec. 2.1). How- 
ever, the LBG magnitude range is more extended towards 
fainter luminosity at all redshifts considered. It follows that 
LAEs represent a luminous LBG subset, whose relative ex- 
tent depends only on the adopted selection criteria. 

• In addition to the previous point, for AIuv ^ —18, —19 
at 2 « 6, 7/8 respectively, LBGs and LAEs share very similar 
values of M* , f * and Aldust , in any given magnitude bin. This 
further strengthens the case for the argument that LAEs and 
LBGs are essentially the same class of galaxies; the differ- 
ences inferred between these two populations solely lie in 
the different techniques employed for their identification. 

• In broad agreement with the results of Stark et al. 
(2010), the faintest LBGs most often show Lya emission 
if and only if the line EW selection criterion is made more 
stringent, i.e. EW > 55A. If instead the threshold EW > 
20A, is used as per the canonical criterion used for identi- 
fying LAEs, we find that all LBGs within the magnitude 
range Muv ~ —23 to —20 would be identified as LAEs at 
2 « 6. Further, the UV range (fraction of LBGs showing 
Lya emission) narrows (decreases) with increasing redshift. 

In spite of the remarkable predictive power of the the- 
oretical model presented, considerable room for improve- 
ments remains. In fact a number of delicate, albeit reason- 
able, assumptions have been made as explained through- 
out the text, which are now summarized. First, the relative 
escape fraction of Lya to continuum photons, fa/fc, used 
is essentially constant at a given redshift with a tiny halo 
mass dependence, that has been fixed by matching the am- 
plitude of theoretical LAE Lya LFs to observations. Such 
ratio sensitively depends on the dust distribution (homoge- 
neous/clumpy) in the ISM of a galaxy. While all existing 
extinction curves (SN, Calzetti, Milky Way) assume dust 
to be homogeneously distributed to find fa/fc < 1, the ra- 
tio can be boosted to fa/fc > 1 if dust is clumped in the 
ISM, as shown by Neufeld (1991) and Hansen & Oh (2006). 
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Modeling of this process is beyond the capabilities of the 
present study since the large cosmological volumes required 
to achieve sufficient galaxy population statistics do not allow 
us to treat fine ISM details. 

Secondly, the observed Lya luminosity also depends 
on the peculiar velocities which have not been included 
in our calculation; inflows/outflows into/from the galaxy 
can blueshift /redshift the Lya line, thereby decreas- 
ing/increasing Ta (Verhamme et al. 2006; Dayal, Maselli & 
Ferrara 2011). However, the extent to which peculiar veloc- 
ities influence Ta is debatable since these calculations have 
mostly been performed under idealized situations. For ex- 
ample, Verhamme et al. (2006) have used spherically sym- 
metric outflows of If I to show an enhancement in Ta ; how- 
ever, many studies, e.g. Fangano et al. (2007) and refer- 
ences therein, have shown that Kelvin-Helmholtz instabili- 
ties would result in breaking-up such symmetric outflows. 

Thirdly, both observationally and theoretically, the de- 
pendence of /esc on the halo mass and redshift have been 
hugely debated with the values obtained by different groups 
ranging between /esc = 0.01 — 0.8 (e.g. Ricotti & ShuU 2003; 
Gnedin et al. 2008; Wise & Gen 2009; Razumov & Larsen 
2010; Vanzella et al. 2010; Nestor et al. 2011). We have used 
a constant value of /esc = 0.02 for all galaxies in each of the 
simulation snapshots, following the results of Gnedin et al. 
(2008) who flnd /esc = 0.01 — 0.03 for galaxies with mass 
^ lO^^M© at z ~ 3 — 5. We remind the reader that /esc 
is a crucial parameter that enters into our calculations, and 
an increase (decrease) in its value would lead to a corre- 
sponding decrease (increase) in the nebular recombination 
line flux, thereby decreasing (increasing) both L^"* and U"^ 
while increasing (decreasing) the size of the Stromgren re- 
gion built by any galaxy. 

Fourthly, due to simulation constraints, we have used 
the stellar mass weighted age and assumed the SFR to be 
constant over such time period to compute the intrinsic spec- 
trum of each galaxy. However, using the complete star for- 
mation history of each galaxy, i.e. using the age, metallicity 
and mass of each individual star particle in each galaxy to 
build the composite spectrum, might affect the results pre- 
sented here, which is the focus of an ongoing study. 

All of these above-mentioned physical effects affect the 
EW calculations presented in this paper; the EW compari- 
son with observations (especially with those of Stark et al. 
2010 carried out in this work) must therefore be seen only 
as very tentative and preliminary. It is hoped that the in- 
creasing complexity and power of numerical simulations, and 
upcoming instruments such as the JWST, ALMA and ongo- 
ing surveys with the HST/Keck will help disentangling such 
issues. 
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